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^ ' Abstract 
m . 

■ Based on the low-energy effective Hamiltonian, we calculate the new physics 
corrections to the branching ratio and the differential distributions of the rare 
decay B ^s77 induced by the new gluonic and electroweak charged-Higgs 
penguin diagrams in the general two-Higgs-doublet model with the restriction 

Q^' Xfj^ = for i 7^ j. Within the considered parameter space, we see the 

^ ■ following: (a) the standard model predictions of B{B — > Xsjj) and ApB have 

a moderate m<j dependence; (b) in model III , the prediction of the branching 
ratio B{B — > ^^77) ranges from one third to three times of the standard 
^ ' model prediction, but is highly correlated with that of B(B Xsj); (c) the 

' new physics enhancement to the branching ratio B{B ^sll) in model II 

can be as large as (30 — 50)%; (d) the contribution from IPR diagrams is 
dominant and hence four normalized differential distributions are insensitive 
to the variation of scale /x and possible new physics corrections; (e) due to the 
smallness of its decay rate and the long-distance background, the B — > Xs^^ 
decay is not a better process in detecting new physics than the B ^ Xs^ 
decay. 



PACS numbers: 13.25.Hw, 12.15.Ji, 12.40.Bx, 12.60.Fr 
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I. INTRODUCTION 



The flavor changing neutral current (FCNC) inclusive transition B X^'-f has been 
a subject of great interest during the past few years. The basic theoretical framework 
is the standard model (SM) at scales or m^. QCD short distance corrections are 
then incorporated via the renormalization group technique to yield the low energy effective 
Hamiltonian valid at scale 0(mf,) which is relevant for B decay processes 

As is well known, the rare decay B — >■ X^'y is theoretically very clear, measured at CLEO 
with increasing accuracy and in remarkable agreement with theoretical estimates. Great 
progress in both the theoretical calculation |^ and the experimental measurement ^ for 
B —>■ Xs7 decay enable one to constrain the new physics models, such as the two-Higgs- 
doublet model (2HDM) 0, the minimal supersymmetric standard model (MSSM) and 
the Technicolor models ||^. In light of these developments, it is natural to consider other 
inclusive channels which as a whole may separate out contributions from various operators 
in the effective Hamiltonian. 

The rare flavor changing neutral current inclusive transition B —>■ Xg'j'j, which is ex- 
pected to be ~ 10^^ smaller in its branching ratio relative to B ^ ^sl, has attracted new 
attention in view of the planned experiments at KEK and SLAG B-factories and existing 
accelerators, which may test branching ratios as low as 10~^. Just like the decay B — >■ Xg'j, 
the decay B — > Xg'yy is relatively clean after some proper precaution to take out effects due 
to strong resonance, for instance, the effect of rjc at its peak to the two photon spectrum 
||T0| . The studies about these rare decay modes will then provide further opportunities in 
testing the whole technology of weak decays, or better yet in providing some clues of new 
physics. 

There have been some theoretical investigations for the process B —>■ Xs'j'j, which corre- 
sponds at the quark level to the transition b —>■ s'-ff. Theoretical calculations in the SM were 
firstly done on the basis of pure electroweak theory [pl] , |T2| and subsequently improved to 
include the leading order (LO) QGD effects |T3|,[T^. Nearly all the studies of B ^ Xg^yy and 
B — > decays are based on the free decay of b quark and can be justified from the heavy 
quark effective theory (HQET). According to the argument in Ref. the HQET correc- 
tions to these rare decay modes are strongly suppressed by powers of (A/m^)^. Furthermore, 
Ghoudhury and Yao []TB| have shown that even when(A/mc)^ and {A/rricY corrections are 
included, the overall HQET corrections to the free quark results of B ^ Xg^y are only a 
few percent and can be safely neglected. In this paper, we will present our results at free 
quark level. 

Recently, we estimated the new physics contributions to the inclusive decay b qg* — *■ 
qq'q' with q G {d,s} and q' G {u,d,s} in the model HI: the third-type two-Higgs-doublet 



model [|l^,|T3], and found that the new physics contributions can be significant [1^. In 



this paper, we calculate the new physics corrections to the branching ratio and differential 
distributions of the rare decay B Xg'j'j induced by the new gluonic and electroweak 
charged-Higgs penguin diagrams in the model HI. 

The plan of this paper is as follows: in the next section, we give a brief review of the 
process B —>■ Xg^j in the SM at LO approximation, present the relevant formulae and 
discuss some general characteristics of this process. In Sec. HI, we will investigate in detail 
this rare decay in the model HI and present some numerical results. In the last section. 
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concluding remarks are added. 



II. B Xsj-f IN THE SM 

By the use of an extension of Low's low energy theorem p!3| , |20[] , or alternatively, by 



applying the equation of motions |21], the most general effective Hamiltonian that describes 



radiative b ^ s decays with up to three emitted gluons or photons is given by Q 

Heff = -^V,,V*J:C^i^^)0^, (l) 

i=l 

where the current-current, penguin, electroweak- and chromo-magnetic dipole operators are 
of the form 







02 = 


{SaJ^LCa){Cf3-f^Lbi3) , 


Os = 


q=u---b 


04 = 


q=u---b 


05 = 


q=u---b 


06 = 


(sa-f^Lbp) J2 (^/37m^?«) , 

q=u---b 


07 = 




08 = 





(2) 



where a and [3 are color indices, a =1, ■ ■ -,8 labels SU{?>)c generators, and L, R={1 =f 75)72, 
while F^y and denote the QED and QCD field strength tensors, respectively. 

The Wilson coefficients Ci{fi) in Eq.(0) are process independent and their renormalization 
is determined only by the basis of operators Oj. In our calculation of the leading order matrix 
elements of 6 — > 577 shown in Fig.|l], we find no regulation scheme dependence enters through 
the new class of penguin diagrams with two external photon. This means that it is suffice 



to use LO regularization-scheme-independent Wilson coefficients Oi,...,6, 07^ and Csg ||22 
and we need not consider the matrix elements due to the insertion of O5 and Oq into the 
one-photon penguin diagrams. In Table |, we present the numerical values of the Wilson 
coefficients at the scale fj, = 0{mb) 0- 



^Strictly speaking, operators {sa'y^Lua){up^fj_Lbp) and {sa^^^Lup){up^fj_Lhci) should be added. 
Here we neglect them for their small coefficients. 

^ These Wilson coefficients are numerically well consistent with those given in Ref. |14|. 
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At the leading order, the amphtude for the decay b{p) — > s{p') + 7(^1) + 7(^2) can be 



expressed as |114 



A 



ie^Gf 



(3) 



with 



F5 

wr 



[N,Ci{fi) + C2{fi)]Qlt^c 

+C3{fi) {n, [Qli^a + Hs + '^b) + Qlit^u + Kc)] + Qlif^s + /tfe)} 

+C4(/i) {Qrf {Kd + Ks + Kb) + QI{Ku + Kc) + NcQI{Ks + Kb)^ 
- [N.C^ifi) + CM] [QliKd + Ks + Kb) + QI{Ku + K, 

[C,{fi) + N,CM]Ql, 
C7iiJ,)Qd , 



m 
1 
2 



■ie^.uxMk2l^ + kl ■ k2g^'' - k'ik^) (1 - 2k,) + 4 U'^'^ - 





(4) 




(5) 




(6) 




(7) 






k,-k2r'' 


(8) 



2p ■ k2 
1 



hl^inibR + msL){^- 1^2 + rJib)-/" 



2p' ■ h 



(9) 



where k\ and k2 are four-momenta of the emitted photons, A^^c = 3 denotes the number of 
colors, Qu = 2/3 and Qd = —1/3 are the up-type and down-type quark electric charges, and 



the factor Kg is defined as 



Kf, 



1 1 (ix , 

- H / — \n[l- ZgX + 

2 Za Jo X ^ 



2 2 



2 ^ z„ 



arctan 



i-Za 



-f + 2 (In 



Za-A 



— 2m In 



/ ^+^Jzq~i 

\ 2 



if Zg<A 



otherwise . 



(10) 



Here Zq = 2ki ■ k2/m?q- The masses appearing in W!/'^ arise from operator O7 while those 
in Wl^^q are internal quark masses in the loops (see Fig.l). In our numerical analysis, we 
take all these quark masses as current mass. The factor k, in Eq.(|l^) is the loop integral 
function p3| and its absorption part reflects intermediate quark threshold effect. In getting 
the expressions of W^'^ that are essential to obtain the matrix elements of 6 — > S77 decay 
presented in |T^, we used the following identities: 

7m7a7<7 = QaXl^L + g^l\'la - g^,alx + i(^i,paxYlb , (11) 

1 Pll fJ-'il fJ-i ~ '^^/il/i2M3A'475 ~ 9p-J,Pi9piP2 ~^ 9n2IJ-i9piHS, — 9h2P--a9piPH "I" 9 p-iPil Pll P2 

~ 9^.2^41^11^3 ~^ 9p2P3.1 fJ-ll fJ-i ~^ 9nipLAlp2lp3 "~ 9 pipsl P2I PA ~^ 9 pi P2I Psi 1^-4 ' (l^) 
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The square amplitude summed over spins and polarizations is then given by 



\A\ 



if 

+2Re[F^F;{l - 2ft:,)]A^5 + 2Re[F^F;]A^^ + 2Re[F^F;{l - 2k,)]A\^] 



rnt{\F2\''A22 + l^spAs + iFyr^r + 2Re[FjF;]A27 



(13) 



where the quantities Aij denote the contractions between the tensor and , and 

the first two terms and the third term of \A\'^ arise from one-particle irreducible (IPI) 
diagrams and one-particle reducible (IPR) diagrams, respectively. In order to give the 
explicit expressions of A^j, we introduced the following notations, 



2ki ■ k2 



2p-k2 



m 



2 ' 



12 



2 ' 



u 



2p-ki 



m 



2 ' 



rhoi = — 7T 



m 



2 ' 



p=l — u — t + s 



P 



m 



2 ' 



U 



(14) 



where rris is s quark mass appearing in Wc^j and p' is the momentum of the outgoing s 
quark. In this framework, coefficients A^j are then given by 



A22 
A55 



A. 



b,s 
25 



A 



27 



A 



77 



/I7 



77 



^77 



^77 



2[(l-p)2-(l + p)s] , 

1 16\Kb\^ + 1(1 - '2Kb)s + 4Kb\^ + Pi[16\ks 



(15) 



[1~2ks)— + 4ks 



s + p) 



+ 16Re {SpiKbK* + s[K,b - 2(1 + Pl)KbK* + PlK*]} , 

±s(l -p^s) , 

on , ^ , (p + pi)g^ + 2(pi -p)tg (p + pi)g2 + 2(pi -p)Mg 
2(1 + pijs H ^ h 



t{s — t) ' u{s — u) 

A57 = Re {%{Kb + PlKs)s - [Api{Kb + Ks) + S [(1 - 2Ks) + Pl(l - 2Kb)]] 



+ 



t(s — t) u{s — m) 

(1 + Pi) [{I- p)A]j + A^ + A 

1 r 2uiu - 2)- 2m 
T 1 + ^+ n ^ 

t 1 — M 1 — M 



(«:;, + piK^)(p-pi)( 



+ 



t(s — t) u{s — u) ' 



3 

77 1 



^ u] 



1-u 



+(P1-P) 



2s 



;i-^)2 

su 



u 



t{l-U) (1-m)2 (1-m)2 



+ (t ^ m) 



-2^t<^ (l + pi)(2 + ut) + 



Pi 



1-p 



2(1 + Pi] 



tu 



m-u) 



2s . . -(1 + pi)(2s + s) - pi(s2 - 2 t u) + 2(2 - p)(l + pi)t m 
' (l-n)(l-t) 



(16) 
(17) 

(18) 



(19) 
(20) 
(21) 



(22) 



(23) 



Unlike Ref. |T^, we here distinguish different s quark contributions which enter in the 
amplitude via p and pi respectively. This is crucial to our results, as described below. Note 
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that in the Ajj term of the square amphtude, there exists infrared (IR) divergence when 
one integrates the amphtude over the physical phase space. This divergence can be canceled 
out when one considers 0{ae) virtual corrections to the b ^ s'y amplitude In order to 



calculate the physical rate of interest that is free of divergence, we have to impose a cut on 
the energy of each photon, which will naturally correspond to the experimental cut imposed 
on the minimum energy of detectable photons. 

In numerical calculations, we use the input parameters listed in Table ^ and take the 
following cuts: 

> 100 MeV , Es > 600 MeV , > , (24) 

9 

where E^ and Eg denote the energy of photon and that of of the outgoing mesons respectively, 
and 6 is the angle between any two outgoing particles. The first constraint is required to 
avoid IR divergence while the last constraint is to exclude photons that are emitted too 
close to each other or to the outgoing s quark. For the mass of s quark, cares must be 
taken. In principle, constituent mass should be used in the phase space integration, while 
the masses appeared in Wj'^ and Wi^'^ should be the current mass. This means that p and 
pi appeared in Eqs. (|T5| - p3D should be 

ml {constituent) rri]^ ml [current) . 

P = 2 - 2 ' Pi = 2 ■ (25) 

m^ mi; 



Following Ref. p5|, the branching ratio of 5 — > Xs'j'j decay can be written as 

th 



r(6 ^ S77) 



X B{B XM)^''^' . (26) 



Using the input parameters presented in Table y, setting p = mf, and taking the cuts given 
in Eq.(l^), we find numerically that 

BiB ^ X,77) ^ 4.6 X 10"^ , ApsiB X,77) = 0.79 , (27) 

where ApB the forward-backward asymmetry of the B ^^77 decay. 

If we use a common s quark mass = 0.5 GeV in the numerical calculations as Ref. 



Ij] did, we find that 

BiB ^ X,77) = 3.8 X 10^^ , ApsiB ^ ^,77) = 0.76 , (28) 



which agree well with the results presented in Ref. fl^, where B{B — > Xs77) = 3.7 x 10 ^ 



and ApB = 0.78. For = 0.15 GeV, however, we find that 

BiB ^ X,77) = 5.1 X 10-^ , ApBiB ^ X,77) = 0.81 . (29) 

From above numerical results, we find that the m^-dependence is weak for ApB, but relatively 
strong for the branching ratio. 

Among all the contributions to the decay rate, the contribution from IPR diagrams is 
predominant, larger than 97% of the total, which is due to the cancellation between Ci and 
C2 and the QCD enhancement of IC7I. Numerical results also show that the branching ratio 
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is sensitive to the cuts we have imposed. For example, if we demand 9 in Eq.(p4D larger 
than 7r/6, the branching ratio will reduce to 3.6 x 10"''. Moreover, because of the scale 
dependence of Wilson coefficients, there exists a ~ 25% theoretical uncertainty at leading 
order approximation, as is the case for B — > Xg'y. 

In Fig.|^, we present the normalized differential distribution {l/r)dr /ds versus s ( solid 
curve ). Comparing with the corresponding result given in [|l^], we find that the peak due to 
cc threshold effect is smeared and the average invariant mass of the two photons is lowered 
when QCD corrections are added. In Figs.(^^ , the normalized differential distribution 
{1 /T)dT /d cos 9^^ versus cos9^^ (sohd curve in Fig.|^) and the spectrum of the two photons, 
defined as the photon with lower energy (solid curve in Fig.^) and the photon with higher 
energy (solid curve in Fig.|^), are plotted. 

From these four figures, one can understand the kinematics of the process as follows: an 
energetic s quark with mean energy around 2.0GeV tends to be emitted, compensated by the 
harder of the two photons; while the less energetic photon tends to go in the direction of the 
s quark. This topology is typical of a bremsstrahlung event of the s quark. Such kinematics 
are very useful for us to separate the short-distance (SD) signal from the long-distance (LD) 



background, which may come from the channel B ^ Xgrj (rj') — > Xs77 |lT0| , p!2| . For example, 
by demanding s < 0.3 and cos 6'^^ < 0, about 75% of the signal is remained while almost all 
the LD background are removed. 

In Figs.(H-^, the contribution from IPI diagrams (dashed curve) and that from the 
interference between IPI and IPR (dot-dashed curve) are also plotted. From these figures, 
we can see that except in some marginal area, the contributions from IPI diagrams and 
the interference are much smaller than those coming from IPR diagrams. Taking this in 
mind, one can infer that though the four normalized differential distributions are calculated 
at leading order approximation, they are insensitive to the variation of scale. The reason is 
that, to a good approximation, a common factor |C7(/i)p can be extracted from both F and 
dr and will disappear in the ratio of them. This feature becomes more evident when the 
cuts of s < 0.3 and cos^^^ < are imposed, which removes nearly all the contributions of 
Oi ~ Oq (see Fig.^]). In fact, a more general statement is that: as long as the new physics 
effects beyond the SM appear only in the matching of the Wilson coefficients of the standard 
effective operator basis, which is so in the MSSM ||25|, and \Cj\ is not suppressed as required 
by 6 ^ S7, these four observable distributions are insensitive to the variation of scale /i and 
the new physics. 



III. B IN MODEL III 



A. Model description 



The two-Higgs-doublet model (2HDM) [0] is the simplest extension of the SM . During 
past years, the models I and II have been studied extensively in literature and tested ex- 
perimentally and the model II has been very popular since it is the building block of the 
MSSM. In Ref. |]12[, the authors studied the B Xg'j'j decay in the SM and models I and 



II in the basis of pure electroweak theory |TT[], and found that the branching ratios in models 
I and II can be appreciably different from that in the SM |T2|. In this paper, we focus on 
estimating the new physics effects on the B ^s77 decay in the framework of the third 
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type of 2HDM, usually known as the model III [|T8|. In the model III, no discrete symmetry 
is imposed and both up- and down-type quarks may couple with either of the two Higgs 
doublets. The Yukawa coupling in quark sector in this case is [0 



(30) 



where 0i and 02 are the two Higgs doublets and rifj^ and are the Yukawa couplings. 
As described in |T^, in order to let 0i correspond to generate fermion masses while 02 
to introduce new interactions, one can choose the following parameterization of the Higgs 
doublets 



'■U,D 



1 

V2 







_]_( V2H+ 
v/2 V + 



(31) 



where v = {V2Gi 



v-l/2 



2A6GeV. After the rotation that diagonalizes the mass matrices 



of quark fields and that of the Higgs doublets, we can get the Lagrangian that is relevant 



for our following discussions. The main features of model HI are as follows [18 



FCNC may exist at the tree level. The neutral and the charged flavor changing cou- 
plings are related by: 



u 

charge 



neutral^CKM , 



'ich 



charge 



^CKMineutral i 



(32) 



where Vckm is the ordinary mixing matrix between down-type quarks ^6], and 



pU,D 
^neutral 



U,D^-I^U,DyU,D 



R 



(33) 



• Like the models I and II, there are also five Higgs bosons in model HI: the light 
and heavy CP-even neutral Higgs boson h'^ and H^, one CP-odd neutral Higgs boson 
and a pair of charged Higgs bosons if^. The transformation relation between 
{H°, H^, H^) in Eq.(|31|) and the mass eigenstates (i/°, A°) can be found in [g. 

In the following of this paper, we will parameterize i^'^trai Ref. ||I8[ 

{^neutral) ij ~ \j ~ 5 (^4) 

and treat as basic free parameters. 



B. Experimental constraints 

There is a considerable interest in the constraint of the parameter space of the 2HDM, 
especially in model HI, since the FCNC may appear at the tree level. Compared with the 
SM, the additional free parameters of model HI are the masses of the additional Higgs bosons 
and the coupling constants Xij. In this subsection, we summarize the main constraints on 
these parameters from direct searches at LEP experiments , and from the measurements 
of F° - FO mixing with F° = S^. 
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Let's firstly turn to mass constraints. The LEP Working Group for Higgs boson searcfies 
28| lias recently reported excess of events that might be due to the production of a neutral 



Higgs boson weighing about 115 GeV, and placed a lower mass limit for the SM Higgs boson 
Mho > 113.5 GeV at 95%C.L.. For the /i° and A° Higgs bosons of the MSSM, the new 
95%C.L. limit ^ is M^o > 88.3 GeV and M^o > 88.4 GeV respectively. But one should 
note that these constraints on the masses of and are not applicable in the 2HDM's 
because the coupling of to Z^Z^ and to A^Z^ go like sin a and cosa( a is a free parameter 
of the model), respectively. For a very small mixing angle a and with A^ [h^) sufficiently 
heavy no lower limit on M^o (M^o) can be set from the LEP data. For the charged Higgs 



bosons in the 2HDM's, the 95%C.L. limit is Mh± > 78.6 GeV ^|. 

Indirect constraints on Higgs masses come from Rb and p, which have been measured 
at LEP. According to studies in Ref. ||3^, if one requires rI^'^°^^ to be within la deviation 
from Rl^^ and p within 2a deviation from p^^^, the preferred range of Higgs boson masses 
in model HI are SOGeV < M^o ~ M^o < UOGeV and Mh± ~ (180 - 220)GeV [0. On the 
other hand, if one allows R^^'^°'^y to vary within 2a errors of Rl^^, a charged Higgs mass of 
a few hundreds GeV is feasible. In addition, unitary condition requires these masses to be 

a lot of processes have to be analyzed to give reasonable 
In Table |T|, we list current constraints on 



18,30 



less than ITeV 

As for the couplings A^^ 
constraints, as have been done in Refs 
these couplings along with the processes from which constraints have been or will be placed. 
It should be noted that despite current constraints on flavor changing couplings A£ and 
are rather loose, some much more stringent constraints will be imposed in the future 
experiments. 

For simplicity, we set (a) all the FC couplings A[^'^(i ^ j) = 0; (b) A^ = Ai for i = u,c, t, 
and = A2 for i = d,s, b; and (c) Ai^2 are real numbers. The advantage of such a setting 
is the new contributions to Wilson coefficients discussed below come only from the charged 
Higgs penguin diagrams and the neutral Higgs bosons are irrelevant to our discussion. The 
model HI therefore differs from the model II only in the couplings of the charged Higgs 
boson to fermions. The Feynman rules for the UiH^Dj and jH^H' vertex are 



—I r 



VcKM ■ i?Al + 75) - ■ Vckm{1 - 75) 



-yH^H : ie{p2 - pi] 



(35) 
(36) 



For more details about the couplings between the Higgs bosons and quarks or gauge bosons, 
one can see Ref. . 

Another very important source of constraint is the process B — > X^'y. In the SM, the 
theoretical prediction of B{B Xg'y) is (2.8 ± 0.8) x 10~^ at LO approximation [^, and 
(3.28±0.33) X 10"^ at NLO approximation The 95%C.L. limit from CLEO measurement 
is 2 X 10""^ < B{B Xs'y) < 4.5 x lO""'. One can see that the theory and the data are 
in remarkable agreement, which implies that only a narrow parameter space of new physics 
can survive. Since the aim of this paper is to estimate the new physics effect on the process 
B Xs'j'j, we will vary the input parameters of 2HDM in a rather large range allowed by 



2.0 X 10-^ BfPiiB X,7) 



< 



2.8 X 10-4 B^ZiB ^ X,7) 



< 



4.5 X 10-4 
2.8 X 10-4 



(37) 
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C. Calculation of B ^ ^s77 

In model III, the operator bases presented in Eq.(^ are insufficient for the process 
B — > Xs77 and we need following additional operators when all the FC couplings are zero 
H: 



0/ = {s^'^^'Rcp){cp-i^Rh^) , 

q=u---b 



OA 



q=u-- 



O' 



167r2 
9s 

IGtt 



(38) 



For the evaluation of Wilson coefficients, we need their initial values with standard 
matching computations. Denoting the Wilson coefficients in the SM with Cf*^(mvK) and 
those from the additional charged Higgs contribution (see Fig.||) with (mw), we have the 
initial values of the Wilson coefficients for the ffist set of operators [ Eq. (EI) ] 



(39) 



where 



Ci'..,{mw) = , 



m 



m 



>^l9i{y) + AiA2^2(2/) + -^>^lgiiy) 



and the explicit expressions of C^'^{mw) can be found in the literature 
Wilson coefficients we have 



C:^^^(m^) = C:^(m^) 



(40) 

For the primed 
(41) 



and C^^ are given as 



Cl^'.imw) = 



= 0, 






1 


2 ^2/1(1/) + 




~ 2 




nib 


1 


X29i{y) + 




~ 2 


"h 


rUb 



m 

rrib 



m 



■5 \2, 



TUb 



(42) 



^Here we neglect the operators Og^io and O9 10 introduced in |35] because in our case, the corre- 
sponding Wilson coefficients keep zero at any scale. 
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In our calculations we neglect the contributions of the internal u and c quarks and define 
functions /ad/), gi{y), and 5(2 (y) as 

= 72(y - 1)3 + T%3T)i ' 

. / X l/(%-3) |/(-3?/ + 2) 
^^(^^=12(^^+ 6(,-l)3 

= 2A{y - 1)3 + i(^3Tji ' 

where y = m1/m\±. Numerical analysis show that fi{y) and f2{y) have the following 
properties: 

/i(y)<0, /2(l/)>0, |/2(i/)|>3|/i(2/),| (44) 

and (/id/) have similar properties. 

Using the renormalization group, we can get the LO Wilson coefficients at any scale. The 
Wilson coefficients C/^^(/i) can be obtained from |^ with Cf^\mw) replaced by Cl^\mw), 
while Cl^^\iJi) are given byQ 

c["i{^^) = , 

C'-j''\l^) = v^C',''\rnw) + ^(r/M - r^'i)C','^\mw) , (45) 

where rj = as(mvi/)/as(At)- 

The square amplitude summed over spin and polarizations can be written as 

\A\jjj=\A\'^ + \A'\'^ + 2\A*A'\ , (46) 

where the first and second term comes from Oi ~ Og and O'l ~ Og respectively and the 
third term is the interference between them. The explicit expression of is the same as 
that in Eq.(|T3p with Ci being replaced by C/^^, and and are given by 



= 7 i^yttv;] mf-^ mF^\A',, + \f;f;\a',, + \f;f^\a',,} , (47) 



^Because the strong interaction preserve chirality, the first set of operators (Eq.(^)) cannot mix 
with its chirality flipped counterparts, the primed operators, and the anomalous dimension matrices 
of the two separate set of operators are the same and do not overlap |3^. This means the evaluation 
of CKn) are the same as in the SM. 
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where 



A' 
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-2 * s - \, ^ + s * (pi - p) * ( ) + {t^u) 



t{s-t) 



t s-t' 



= Re 



+4:{Hb + Ks)(pi - P) 



U 



S-t 



r. ^ .1 .9 1 / 2 2 2(1 -pi - s) 
A'^ = 2 * (1 - p)AIj + A^, - 2s 2 + - + - + ^ 

(pi -p)(3 + 5pi) - (l + pi)g _ (l-pi)(s^ + 4pitM) ^ 

(s — t){s — u)tu 



(s — )f:)(s — u) 



(48) 



In Eg ■ (|47| - ^8D , is the current mass of s quark and the exphcit expressions of Ajt,A\-j 
and A^^ have been given in Eqs. (|20| - p2D . For l^lf//, two characters should be noted. Firstly, 
the interference term between Oj and 0[ is suppressed by nis/mb and therefore, to a good 
approximation, can be neglected. Secondly, the total IPR contributions are proportional 
to (I-F7P + \F^\'^)Atj. Since \F^\'^ + [F^p is not suppressed, as required hy B ^ in 



2HDM [^, the IPR contribution is then predominant, which is very similar with the case 
in the SM. As we discussed in section II, these characters mean that the four differential 
distributions defined in last section cannot have a large deviation from SM predictions and 
we do not expect to see any new physics signal by the measurement of these four differential 
distributions. Indeed, we have checked our numerical results in a vast parameter space 
constrained by Eq.(0), and find it correct. 

We are now ready to calculate the new physics corrections to the decay rate of i? — >■ Xg^^. 
In order to reduce the effects of theoretical and experimental uncertainties, we define two 
ratios 



R-y-tip) 



»5M 



(49) 
(50) 



where the four branching ratios denote the LO theoretical predictions in the SM and model 
III, respectively. 

In the numerical calculation, we consider the following three typical scenarios. 

• Scenario I: |Ai| = IA2I = A and m^^i is about a few hundred GeVs. This scenario 
can be further divided into two limiting cases: 



Case-1: Ai 
Case- 2: Ai 



A2 = A. 

— A2 = A or A2 



-Ai = A, where A is positive. 
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• Scenario II: A2 = —2 / Ai = A and mH± is about a few liundred GeVs. By setting 
A = -^/2 tan/?, one reproduces tlie numerical results of the Model II, a very popular 
version of the two-Higgs-doublet models. 

• Scenario III: IA2I ^ |Ai| and is about 200 GeV, which is favored by current 
experiments. 

Let us firstly concentrate on Scenario I. From the explicit calculations, we find that the 
new physics contribution in case-1 tends to increase the value of Cj. In Figs. and |^, we 
show the A dependence of the ratio R^^ in case-1 for fixed mH± = SOOGeV and assuming 
/i = 2mfe (dashed curve), nib (solid curve) , mb/2 (dotted curve), respectively. In determining 
the range of A, we have put the constraint of Eq.(p7D and required C7 to be negative and 
positive respectively. 

For the case-1, one can see from Figs. |^ and |^ that 

• The /i dependence of the ratio R^^ is weak for a negative C7: we have 0.7 < R^^ < 1.0 
within the range of 1/2 < fi/rrii, < 2. 

• The /i dependence of the ratio R^^ is strong for a positive C7: we have 0.3 < R^/^/ < 2.7 
within the range of 1/2 < n/rrib < 2; which implies that the new physics contribution 
to the branching ratio B Xg'j'j can be significant. 

• In the same region of A as specified in Figs. ^ and |^, the numerical results show that 
the ratio R^y/Ry only varies within a very small range of 1 ±0.02 for 1/2 < fi/mi, < 2. 
This fact reflects the dominance of IPR contribution to the decay rate of B —>■ Xg'j'y, 
and a strong correlation between the branching ratios of B ^ ^s77 and B Xs'y. 

The main feature of case-2 is that the new physics contribution tends to decrease the 
value of Cj, as a result, coefficient C7 keeps negative for all the values of A. In Fig. ||, we 
show the dependence of R^^ on A in case-2 for fixed mH± = 300GeV and /i = 2mb (dashed 
curve) , nib (solid curve), nib/2 (dotted curve). It is easy to see that the dependence of 
R^-y is weak. Study of the ratio R-y-y/R^ gives the same conclusion as in case-1. 

In Fig. 1^, we plot the niH± dependence of the ratio R^^y in Scenario I. As one can see 
from this figure, R^^ tends to approach 1 when the charged-Higgs boson becomes heavier, 
which is just the decouphng behaviour of the new heavy particle. 

Secondly, we consider the Scenario II. By setting A = — v^tan/5, one reproduces the 
result of the popular model II: the second type of the two-Higgs-doublet models. Wilson 
coefficient Cj keeps negative for all the values of A. In model II, two additional free parame- 
ters will enter into our calculation: the charged Higgs mass nif{± and the ratio tan/3 = V2/V1 
with V2 and Vi are the vacuum expectation value of the two Higgs doublets 0i and 02- 



In Fig. 11, we plot the tan/3 and /i dependences of the ratio R^^ in model II for fixed 
nijj± = SOOGeV^. In this figure, the three curves correspond to /i = 2nib (dashed curve), nib 
(solid curve), and mb/2 (dotted curve), respectively. One can see that (a) the ratio R^^ is 



-"The A defined in this paper is different from that in Ref. IsO] by a factor of \/2- 
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larger than 1.32 in the whole considered range of tan/3 and fi, and (b) the ratio R.y^ has a 
weak tan/? dependence for tan/3 > 3, but sensitive to the variation of scale fi. 

In Fig.|T2|, we plot the tan f3 and mH± dependences of the ratio R^^ in model II for fixed 
fi = rrih. In this figure, the three curves correspond to tan/? = 1 (solid curve), 10 (dashed 
curve), and 100 (dotted curve), respectively. One can see that the ratio R^-y is insensitive 
to the variation of tan (3. 

Finally, we consider the Scenario III. The main characteristic of this Scenario is that C7 



increases with increasing Ai for fixed A2. In Figs. |T3| and |1J, we illustrate the Ai and fi 
dependences of the ratio R^.y for fixed mH± = 200 GeV and A2 = 50. In these two figures, 
the three curves correspond to /i = 2mfe (dashed curve), nib (solid curve), and nib/ 2 (dotted 
curve), respectively. In determining the range of Ai, we have considered the constraint of 
Eq. ( p7D and required C7 to be negative and positive for Figs. |1^ and |T^, respectively. Again, 
the ratio R^^/R^ only varies within a very small range of 1 ± 0.02 for 1/2 < n/mb < 2. 

In Fig. |T3|, we show the A2 and /x dependences of the ratio R-y^ for fixed mjy± = 200 
GeV and Ai = 0.2. In this figure, the three curves correspond to /i = 2mb(dashed curve), 
nib (solid curve), and nib/ 2 (dotted curve), respectively. In determining the range of A2, 
we have required Eq.(^) to be satisfied. It is easy to see that the ratio R^^ has a strong 
(moderate) dependence on fi (A2). The future NLO calculation of R^^ in the SM and 2HDM's 
will decrease the fi dependence. Study of the ratio R^^/R^ also shows strong correlations 
between the branching ratios of B ^ Xsj and B ^^77 decays. 

In Figs. and |T^, we plot niH±, Xi and A2 dependences of the ratio R^^ in model III 



and assuming /i = nib. In Fig. |16|, the soild and dashed curve refers to Ai = —0.018 and 
0.010, respectively. When choosing the value of Ai, we keep C7 to be negative and require 
IC7I > IC^'^^ ( sohd curve) and IC7I < |Cf*^| (dashed curve), respectively. One can find 
from Fig. that the ratio R^^ tends to approach 1 when charged Higgs boson is becoming 
heavier. In Fig. |l^, the three curves correspond to Ai = 0.2, A2 = 40(dashed curve), 50 
(soild curve), and 60 (dotted curve), respectively. To determine the range of mj^-i, we have 
required Eq.(|37|) to be satisfied and C7 > 0. Numerical results also show that for a much 
larger niH±, Cj is driven to be negative. 



IV. CONCLUSION 

In this paper, based on the low-energy effective Hamiltonian, we calculated the rare decay 
B — s> Xs77 in the SM and the general two-Higgs-doublet model with the restriction }^-^ = 
for i ^ j. We focused on the estimation of the new physics contributions to the branching 
ratio B{B —>■ Xg'yy) and the differential distributions {l/r)dr/ds, {l/r)dr/dcos6^^ and 
{l/T)dT/dX^. 

Within the considered parameter space allowed by currently available data, we found the 
following: 

• In model III, the prediction of the branching ratio B{B ^s77) ranges from one third 
to three times of the SM prediction, but it is highly correlated with the corresponding 
theoretical prediction of B{B —* Xg'y). 

• In model II, the new physics enhancement to the branching ratio B{B —>■ ^^77) can 
be as large as (30 — 50)% with respect to the SM prediction. 
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• In the SM and model III, the contribution from IPR diagrams is dominant and hence 
those four observable differential distributions are insensitive to the variation of scale 
11 as well as possible new physics corrections considered in this paper. 

• Although the process B — > ^^77 provides many new physical observables, it is not a 
better process to detect new physics than B — > ^^7 because of the smallness of its 
decay rate and the long-distance background. 

At the next-to-leading order, the coefficients of O3 — Oq and O'l — Og may get enhanced. 
However, due to their subleading feature, we do not expect drastic deviation from the LO 
predictions. 
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TABLES 



TABLE I. Values of the regularization-scheme-independent LO Wilson coefficients Ci(/x) using 
the input parameters as listed in Table II. 





Ci 


C2 


Cs 


C4 


C5 


Ce 


C7 


mb/2 


-0.315 


1.143 


0.015 


-0.032 


0.009 


-0.041 


-0.340 


rrib 


-0.225 


1.10 


0.010 


-0.023 


0.007 


-0.028 


-0.304 


2mb 


-0.154 


1.061 


0.007 


-0.016 


0.005 


-0.019 


-0.273 



TABLE II. Valiu;s of the iiipul parameters uscxl in the; miuierical ealculations. 



as{mz) 




mt 


mz 


mw 


nib 


rUc 


0.118 




llhGeV 


91.2GeV 


80.4Gey 


4.8GeV 


l.bGeV 


ms{current) 


rriK 


rriu 




\VtlVtb\ 


\Vbc\ 





0.15GeV 0.5GeV b.lMeV 9.0MeF 0.04 0.04 0.11 



TABLE III. Current constrains on \ - along with the processes from which constrains have 



been or will be placed. Couplings that do not appear in this table are not constrained. 



coupling 


current constrain 


current constrain from 


future constrain from 




< 1 


— mixing 




^bd 


< 1 


B,i — mixing 




\^ 


< 1 


]jQ _ /)0 uiixing 




\t 


< 0.5 


Rb, Bd - B^i mixing 






> -10 


Rb. R,: 




^sb 


< 40 


b — ^ sec and Bg — Bg 


Z bs, Bl^ 1+1-, b ^ X,n+ii- 








e^e~ tc, t ^ cj, cZ, eg 
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FIGURES 




FIG. 1. Examples of Feynman diagrams that contribute to the matrix elements < s^'-f\Hf,ff\b >. 
The IPI diagrams illustrate possible insertion of Oi ^ Oq, while the IPR diagrams represent the 
insertion of O7. 




_ I I I \ I \ I \ I 

0.0 0.2 0.4 0.6 0.8 

s 

FIG. 2. The normalized distribution l/TdT/ds versus s in the SM (solid line). The contribution 
from IPX diagrams ( dashed line) and that from the interference between IPI and IPR diagrams 
(dot dashed line) are also shown. In plotting this diagram, we set fx = mi,. The mean value of s is 
0.18 for the cuts as given in the text. 
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O 
O 




-- lOOxlPI 
- lOOxInt. 



_ I \ \ \ \ \ \ t \ t 

-1.0 -0.6 -0.2 0.2 0.6 1.0 

FIG. 3. The normalized distribution l/rdr/(i(cos 0^^) versus cos0^^ in the SM (sohd curve) 
for fixed [i = mi,. The contribution from IPI diagrams (dashed curve) and the interference between 
IPI and IPR diagrams (dot-dashed curve) are also plotted. The mean value of cos 9^^ is -0.70 for 
the cuts as given in the text. 

7 I ' 1 ' 1 ' 1 ' 1 ' 1 ' 1 ' 1 ' 1 ' 1 ' 



lOOxlPI 



lOOxInt. 




0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 

FIG. 4. The spectra of lower energy photons (solid curve) in the SM. The contribution from IPI 
diagrams (dashed curve) and that from the interference between IPI and IPR diagrams (dot-dashed 
curve ) are also shown. In plotting this diagram, we set fj, = mi,, X-y = 2E^/mi, and normalize the 
spectra to the total QCD corrected rate. 
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FIG. 5. Same as Fig.4, but for the spectrum of high energy photon. 



1^9 

) / ^ \ 
b } I u,c,t I 



u,c,t I 



/ , \ 

b , U,C,t ^ ) S 



1,9 



u, c, t 




FIG. 6. The Fcynman diagrams relevant to the decays b 57 and b ^ sg m model III. Prom 
these diagrams, C/^ and C^' ^ can be extracted. The internal quarks are the upper type u, c and 
t quarks. 
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FIG. 7. Rryy 3s 3j function of the coupling A in the Case-1 of Scenario I for fixed m/^± = SOOGeV. 
The dashed, sohd and dotted curve correspond to fi = 2mh, ruh and mb/2, respectively. In deter- 
mining the range of A, we have required Eq.(p^) to be satisfied and Cy to be negative. 
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FIG. 8. Same as Fig. 7, but for a positive C7. 
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0.5 I ' ' ' ' ' ' ' ' ' ' ' ' ' ' 1 

250 350 450 550 650 750 850 950 

FIG. 10. The dependence of R^^ on the charged Higgs mass mjj± for Case-1 (sohd curve) and 
Case-2 (dashed curve) in scenario I. In either case, C-j is negative and we set A = 1 and ji = mi,. 
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1 .7 




^ 2 I ' ' ' ' ' ' ' ' ' 

10 20 30 40 50 

ta nP 

FIG. 11. tan/3 and fi dependences of the ratio R-y-y in model II for fixed mff± = 500 GeV. The 
dashed, solid and dotted curves correspond to ^ = 2m5, rub, and mb/2, respectively. In numerical 
calculation, we have required the constraints as given in Eq.(^) to be satisfied. 




^ Q I I \ I \ I \ I \ I \ I \ I 

'300 400 500 600 700 800 900 1000 

rriH 

FIG. 12. mjj± and tan/3 dependences of R^^ in model II for fixed /i = nib- The solid, dashed 
and dotted curves correspond to tan/3 = 1, 10, 100, respectively. The last two curves can not be 
separated clearly. In numerical calculation, we have required the constraints as given in Eq.(p7|) to 
be satisfied. 
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0.01 5 



FIG. 13. The dependence of i?^^ on Ai in scenario III for fixed A2 = 50 and mf{± = 200GeV . 
The dashed, sohd and dotted curves corresponds to ;U = 2m;,, m^and nib/ 2 respectively. In deter- 



mining the range of Ai, we required Eq.(37) to be satisfied and C7 to be negative 




FIG. 14. Same as Fig. 13, but for a positive C7 
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0.0 I ' ' • ' • ' • ' • ' • ' 

40 42 44 46 48 50 52 

A/2 

FIG. 15. The relationship between R^^ and A2 in scenario III for fixed Ai = 0.2 and 
rrifji = 200GeV. The dashed, solid and dotted curves corresponds to /i = 2mb, mf,and mf,/2 
respectively. In determining the range of A2, we have required Eq.(|37[) to be satisfied and A2 3> Ai. 



1 .6 




0,8 - 

Q g I I \ I \ I \ I \ I \ I \ I \ I \ I 

'100 200 300 400 500 600 700 800 900 1000 

FIG. 16. R-y-y as a function of the charged Higgs mass mfj± for A2 = 50 and fi = rrib ■ Solid 
and dashed curve correspond to Ai = —0.018 and 0.010, respectively. When choosing the value of 
Ai, we keep Cj to be negative and require IC7I > jCf^^l (sohd curve) and IC7I < ICf*^! (dashed 
curve), respectively. 
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0.5 




350 



m. 



FIG. 17. mfj± dependence of R-y-y for fixed Ai = 0.2, ^ = rub and A2 = 40 (dashed curve), 50 
(solid curve) , 60 (dotted curve) . To determin the range of mff± , we require Eq. ( p7| ) to be satisfied 
and C7 > 0. Numerical results also show that for a large mjj±, C7 is driven to be negative. 



27 



